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ABSTRACT 


The effects of aerodynamic loads induced by static aeroelasticity 
are analyzed for the Saturn configurations. This aeroelastic load 
growth is shown, for any given Mach number, to increase at an acceler- 
ated rate with increases in dynamic pressure. During the period of 
high dynamic pressure, the Saturn control requirements due to aero- 
elastic effects are significant. The unexplained external body moments 
recorded during the Saturn I Block II flights are c.orrelated to aero- 
elastic load growth. Local aerodynamic bending moments are increased 
at a time in flight when the vehicle loadings are already large and 
near design limits. This aeroelastic load growth is additional to that 
of flexible bending mode dynamics and should be included when analyzing 
wind profiles for vehicle loadings. 
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DEFINITION OF SYMBOLS 


Symbol 


Definition 

AjB j C j I**) 

I 

curve fit angle-of-attack distribution coefficients 
representing the aeroelastic load growth 

3. , j b , , ••• 9 

3 3 

i . 
3 

elastic curve-fit coefficients corresponding to an 
equilibrium angle-of-attack distribution of a = qJ . 

a^. ,b^_, • • . j 

4 

elastic curve-fit coefficients corresponding to 
the static angular deflection of initial loading 

CP/D 


rigid body center of pressure, calibers from sta- 
tion 100 

£CP/D 


aeroelastic load growth center of pressure, 
calibers from station 100 

CG/D 


vehicle center of gravity, calibers from station 
100 

% 


rigid normal force coefficient gradient, per degree 



incremental normal force coefficient gradient due 
to aeroelastic load growth, per degree 

( C N a ) FLEX 


flexible normal force coefficient gradient (sum of 
rigid and aeroelastic), per degree 

Ci 


aerodynamic moment coefficient about c.g., 
(rad/sec 2 ) (1/ rad) 

D 


vehicle reference diameter, meters 


M AEROELASTIC 
N* /M 


aerodynamic moment about vehicle's center of 
gravity due to the aeroelastic load growth, N-M 

aerodynamic force coefficient, (m/sec 2 ) (1/rad) 


q 


free stream dynamic pressure, N/M 2 


S 


vehicle reference area, M 2 




initial angle-of-attack distribution over vehicle, 
radians 


iii 



DEFINITIONS OF SYMBOLS (Continued) 


Symbol 


Definition 


°E 


elastic angle-of-attack distribution over vehicle, 
radians 


°F 


final angle-of-attack distribution over vehicle, 
radians 


T] 


nondimens ional vehicle station, where tj = 0 at gimbal 
station 


Superscripts 

' elastic curve-fit coefficients based on a unit value 

of dynamic pressure, q 
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TECHNICAL MEMORANDUM X-53634 


AEROELASTIC LOAD GROWTH EFFECTS ON SATURN CONFIGURATIONS 

SUMMARY 


The static aeroelastic effect on aerodynamic loadings was deter- 
mined for the Saturn configurations. The results indicate an appre- 
ciable increase in the vehicle bending loads. For the Saturn V 
configuration, vehicle station 3100 inches, a critical stage splice, 
had a 26 percent increase in shears and bending moments due to static 
aeroelasticity at a flight condition of 12 degrees angle-of-attack at 
Mach 1.2. This 26 percent increase is composed of 19 percent due to 
static deflections and 7 percent due to iterative aerodynamic load 
growth. The net aeroelastic load growth varies along the vehicle 
length and, consequently, causes changes in the total center of pres- 
sure which alter the vehicle's stability. 

The nonlinearity of aeroelasticity with aerodynamic lift distri- 
butions and dynamic pressure precludes any readily available numerical 
factor by which this effect can be superficially included. In general, 
as the vehicle flight loads are increased, the aeroelastic effect also 
increases. Saturn flight tests have already indicated many discrep- 
ancies in aerodynamic data. This flight test error, hitherto labeled as 
an unknown external aerodynamic moment, is attributed *:o static 
aeroelasticity. 


I. INTRODUCTION 


The effect of elastic vehicles on flight loads has, in the past, 
been limited to the analysis of bending mode feedback through the closed 
loop control sensors. While such an approach includes vehicle bending 
mode dynamics, there is no provision for determining the aerodynamic 
load growth induced by vehicle flexibility. This report presents an 
analytical evaluation of the aerodynamic load growth or steady state 
aeroelastic effect for Saturn configurations. 

Aeroelastic considerations are a particularly important factor in 
determining the steady state flight loads of a vehicle because, in a 
flight environment, the vehicle's deflected shape causes changes in 
local angle of attack. These changes in local angle of attack will 
induce an additional aerodynamic loading causing further increases in 
local angles of attack. The resultant deflected shape is an equili- 
brium between local aerodynamic forces and vehicle stiffness. 



A noniterative superposition method of aeroelastic analysis is pre- 
sented in this report. Although the analytic evaluation was based on 
quasi-steady aerodynamics, the nonlinear aerodynamic characteristics of 
the lift gradient were included in the analysis. 


II. ANALYSIS 


The steady state aeroelastic effects on flight loads are evaluated 
from the angle-of-attack distribution at which equilibrium is established 
The mathematical approach presented here is based on a noniterative super 
position technique which is now described. 

The following angle-of-attack symbols are used: 

initial angle-of-attack distribution 

elastic angle-of-attack distribution 

final angle-of-attack distribution at which equilibrium is 
established. 

Their relationship is defined as 


Oj. = 
C£p = 


Op CKj + Og. 


( 1 ) 


The problem is to determine the elastic angle-of-attack distribution for 
any given flight condition. For the Saturn vehicle, the elastic angle- 
of-attack distribution can be sufficiently represented by an eighth order 
polynomial of vehicle station. Letting q be the nondimens ional vehicle 
station, where q = 0 at the gimbal station and t) = 1.0 at the nose, then 

Og = A + Bq + Cq 2 + Dq 3 + Eq 4 + Fq 5 + Gr| S + Ht) 7 + Iq 8 . (2) 


To develop this aeroelastic solution, it is necessary to digress at this 
point and contemplate known solutions. For example, if the equilibrium 
angle-of-attack distributions were known, then the elastic angle-of- 
attack distribution could readily be computed as static angular deflec- 
tions from the vehicle's structural properties and the known loadings. 
The results can then be curve-fitted, in this instance, to an eighth 
order q polynomial. Therefore, if 

a F j = V < 3 > 


2 



where j = 0, 1, 8, then, the curve-fitted polynomial coefficients 

can be computed for each elastic angle-of-attack distribution. 

OL. = a + b q + c q 2 + d q 3 + e q 4 + f .q 5 + g .q 6 + h q 7 + i,q 8 . (4) 

EJ J J J J J J J J J 

For the general aeroelastic solution, the initial angle-of-attack 
distribution due to the static deflection can be computed and represented 
as an eighth degree q polynomial: 

= a l + ^1^ + C I^ 2 + ^i 1 ! 3 + e x^ 4 + + + ^ + ^x 1 ! 8 * (5) 

Now substituting equation (2) into equation (1), we obtain 


Cdp = Oj. + A + Bq + Cq 2 + Dq 3 + Eq 4 + Fq 5 + Gq 6 + Hq 7 + Iq 8 (6) 

and, by using the definitions of equation (3), 


Op = a T + Aa p0 + Bo^ + COp 2 + Da^ + EOp 4 + Fo^ + GOp 6 + HOp 7 + IOpg. 

(7) 

Now, writing the initial angle-of-attack distributions for each term in 
equation (7) yields 

Of = Oj. - (a ] . + b x q + Cj.q 2 + d.j.q 3 + e ]; q 4 + f^ 5 + g^ 6 + l^q 7 + i ][ q 8 ) 

+ Ba + ca l2 + d« i3 + Ea l4 + Fa I5 + Ga l6 + Ha I7 + ia Ig , (8) 

which will reduce to 

Aa l0 + Ba^ + cai I2 + . . . + la Ig = a 1 + bj-h + ^q 2 + • • • + i^n 8 * ( 9 ) 
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Now the initial angle-of-attack distribution of each term on the left- 
hand side of equation (9) can be written as the final angle of attack 
minus the elastic angle-of-attack distribution. 


A [ (1 - a Q ) - b Q q - c Q q 2 - . . . i Q q 8 ] + B[-a x + (1 - b x ) q - c x q 2 


- ... - ixq 8 ] + C[-a 2 - b 2 q + (1 - c 2 ) q 2 - ... - i 2 q 8 ] 


+ ... + I[-a a - b a q - c 8 q 2 - . . . + (1 - i a )q 8 ] 

= CKj. + bj.q + c x q 2 + d I q 3 + e ] .q 4 + f ] .q 5 + g^ 6 + h ] .q 7 + g^ 8 . (10) 


By equating power coefficients of q, the coefficients of the elastic 
angle-of-attack distribution can be determined from the following 
equation: 
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(ID 

It can be seen that the basic equation for the aeroelastic coef- 
ficients, equation (11), bears much similarity to the characteristics 
equations for vehicle bending dynamics. There are, in fact, many more 
subtle similarities, including divergent vehicle velocity, which is 
analogous to a natural frequency in bending dynamics. This nonlinearity 
with dynamic pressure is more evident by defining a set of prime coef- 
ficients which represent the deflected curve fits based on a unit value 
of dynamic pressure. 


4 



Then 


a = a' q 
n n 


i 

n 



( 12 ) 


where n = 0,1,2, ..., 8, and I. Now, dividing both sides of equation (11, 
by q results in the following set of characteristic equations: 



(13) 


This technique can be extended to include the aeroelastic fin twist 
effect (see reference 1). 
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III. RESULTS 


A. Aeroelastic Effect on Saturn I Block II 
and Comparison with Flight Tests 

The aeroelastic load growth characteristics of the Saturn I Block II 
vehicles were determined for various Mach number and dynamic pressure com- 
binations. Quasi-steady state aerodynamics was assumed, and the aero- 
dynamic loadings for the various angle-of-attack distributions [2], as 
defined by equation (3), were computed by multiplying the local lift 
gradient with the local angle-of-attack distribution. The coefficients 
of the final angle-of-attack distribution (A -» I, equation (11)) were 
evaluated, and the aerodynamic loading on the vehicle determined using 
superposition techniques. 

Since the equilibrium loadings are referenced to that of a unit 
angle of attack at the gimbal station, the vehicle's deflected shape 
increases the lift coefficient and moves the center of pressure forward 
of that predicted by theoretical aerodynamics for an undeflected vehicle. 
In figure 2, the increase in lift coefficient gradient, ACm , is seen 
to be dependent on both Mach number and dynamic pressure. This is not 
too surprising since the aeroelastic characteristic equations (equation 
(11)) are nonlinear with dynamic pressure. The center-of-pressure loca- 
tion of this lift gradient is indicated as A3P/D in terms of calibers 
forward of the gimbal station (figure 3). The aerodynamic moment about 
the vehicle's center of gravity due to this aeroelastic load growth is 

M . . = (Z% ) - yf) OqSD. (14) 

aeroelastic v w q; \ D Dy 


Flight tests results of Saturn I Block II vehicles confirm the 
analytical calculations of the aeroelastic load growth. The flight 
instrumentation of these vehicles included both Q-Ball and fin-mounted 
angle-of-attack meters. With no aeroelastic load growth considered, 
flight test correlation between the angle-of-attack meters and the 
gimbal engine deflections required an external applied moment. This 
unexplained moment, which was described as aerodynamic in nature, was 
observed to be relatively large in pitch and yaw. The peak values of 
these unexplained external moments, as recorded in the MSFC Flight 
Test Evaluation Reports [3, 4, 5, 6], are shown in the last column of 
Table I. Using the flight conditions at the time of the peak external 
moments, the aerodynamic moments due to the aeroelastic load growth 
were computed from equation (14). Comparisons between the observed 
flight test unexplained external moment and the moment due to aero- 
elasticity indicate good agreement in magnitude and direction for all 
flights in both pitch and yaw. Thus, inclusion of the aeroelastic load 
growth would have greatly reduced or eliminated the unexplained external 
moments from reconstructed flight simulations. 
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TABLE I 


Comparison of Aeroelastic Moment With 
Observed Flight Test External Moment 


Saturn 

Flight 

Vehicle 


Time 

(sec) 


a 

(deg) 


SA-7, Pitch 76.5 
SA-7, Yaw 72.0 
SA-8, Pitch 64.0 
SA- 8 , Yaw 72.7 
SA-9, Pitch 63.5 
SA-9, Yaw 69.0 
SA-10, Pitch 75.0 
SA-10, Yaw 75.6 



^N 

Z£P 

Moment Due 
to Aero- 

Observed' 
Flight Test 

D 

elasticity 

(N-M) 

Unexplained 
External 
Moment (N-M) 

.0134 

5.12 

376,000 

700,000 

.0128 

4.88 

430,000 

430,000 

.0095 

4.53 

268,000 

650,000 

.0115 

5.02 

205,000 

270,000 

.0105 

4.64 

316,000 

650,000 

.0109 

4.77 

349,000 

270,000 

.0115 

5.05 

247,000 

500,000 

.0112 

4.96 

246,000 

200,000 
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Since (£CNq;) aeroelastic *- s a f unction of both Mach number and dynamic 
pressure, the final flexible lift gradient and center-of-pressure coef- 
ficients must be evaluated for each specific flight trajectory. Pos- 
sibly the best scheme would be to incorporate equations (15) and (16) 
within the program and use a look-up table as part of the input for 
(■^No^aeroelastic anc ^ (f£P/D) ae roelastie (figures 2 and 3), 


B. Local Bending Moment Effect on Saturn IB 

The static aeroelastic effect on local aerodynamic bending moments 
was investigated for the Saturn IB configuration (figure 4). The per- 
centage of increase of the induced aeroelastic bending moment component 
to the rigid body aerodynamic bending moment was plotted as a function 
of vehicle station for various Mach numbers with a fin angle of attack 
of one degree and dynamic pressure equal to 20,000 and 30,000 N/M 2 in 
figures 5 and 6, respectively. The aeroelastic effects increase slightly 
from the gimbal station toward the nose. For q = 20,000 N/M 2 , the aero- 
elastic increase is between 7 and 14 percent, while for q = 30,000 N/M 2 , 
the aeroelastic increase is between 11 and 22 percent. 

To determine the aeroelastic nonlinearity with angle of attack, a 
particular vehicle station (station 1663 inches) was analyzed for various 
angles of attack and dynamic pressures at Mach 1.0 and Mach 1.2 [7], 

The results shown in figures 7 and 8 indicate the nonlinearity angle- 
of-attack effect to be significantly more important at Mach 1.0 than 
Mach 1.2. Also, for any given condition, the percentage of local aero- 
dynamic bending moment will increase with dynamic pressure. 


C. Local Bending Moment Effect on Saturn V 

An aeroelastic analysis was made for Saturn V using linear aero- 
dynamics in the low angle-of -attack region (a < 2°). Figure 9 summar- 
izes the percentage of increase of the local aerodynamic bending moment 
at station 3100 inches resulting from aeroelastic load growth. At 
increasing values of dynamic pressure, the aeroelastic component increases 
in significance. These percentages of increase in local aerodynamic bend- 
ing moment of Saturn V are reasonably similar to those of Saturn IB equiva- 
lent vehicle station 1540 inches (figures 5 and 6). 

Using the dynamic pressure profile consistent with the standard 502 
flight trajectory, the aeroelastic load growth was computed using non- 
linear aerodynamics at various angles of attack for Mach numbers 0.9, 

1.0, and 1.2 [8]. The resulting increases in aerodynamic bending 
moments at vehicle station 3100 are shown in figure 10. A load condi- 
tion of 12 degrees angle of attack at Mach 1.2 will induce approximately 
26 percent increase in the aerodynamic bending moment. Nineteen percent 
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of this 26 percent aeroelastic increase in bending moment is due to 
static deflections and 7 percent is due to iterative aerodynamic load 
growth. 

The aeroelastic effects on the Saturn V normal aerodynamic force 
coefficients and aerodynamic moment coefficients were evaluated at 
Mach 1.2 (figures 11 and 12) and at Mach 1.7 (figures 13 and 14). For 
comparisons, these total vehicle coefficients were computed (1) for a 
rigid vehicle (shown as solid lines), (2) for a flexible vehicle bend 
under original loading (shown as dashed lines), and (3) for a flexible 
vehicle whose deflected shape results from considering both bending and 
iterative load growth, or total aeroelastic effect (shown as dotted 
lines). Examination of these total vehicle aerodynamic coefficients 
indicates that the total aeroelastic effect is significant with respect 
to the bending effect particularly at large dynamic pressures and large 
angles of attack. At Mach 1.2 with a dynamic pressure of 30,000 N/M 2 
and a 12-degree angle of attack, bending effects increase rigid body 
aerodynamics by 9 percent in normal force coefficient and 35 percent in 
the total aerodynamic moment coefficient. Analysis of the same condi- 
tion, using the total aeroelastic effect, which includes both bending 
and iterative load growth, results in increases to rigid body aero- 
dynamics amounting to 12 percent in normal force coefficient and 43 per- 
cent in total aerodynamic moment coefficient. Similar comparisons are 
indicated at Mach 1.7 in figures 13 and 14. 

The bending and total aeroelastic evaluations for the aerodynamic 
normal force and moment coefficients shown in figures 11 through 14 
were obtained using an eighth order polynomial curve Irt solution for 
the vehicle's deflected distribution. A similar approach could, in 
fact, be obtained from equation (11) using generalized bending mode 
coordinates. Since these coordinates are a convenient means of intro- 
ducing vehicle bending dynamics into trajectory control studies, their 
use in evaluating aeroelastic characteristics was investigated. The 
aerodynamic coefficients resulting from changes in local angle of attack 
due to bending and static aeroelas ticity were computed using an eighth 
degree polynomial and compared with results using as many as the first 
four bending modes. In addition, an exact solution of these aerodynamic 
coefficients was obtained from actual angular deflection resulting from 
the initial loading and the final equilibrium loading. The results of 
this comparison (see table II) indicate that considerably more than four 
bending modes are required to approach the accuracy of an eighth order 
polynomial. Furthermore, by comparing the bending column of bending 
mode approximations to the total aeroelastic value of the exact solution, 
yet neglecting the iterative aeroelastic load growth, the effect of 
including bending modes can be evaluated. For this particular flight 
condition, including only bending effects as approximated by the first 
bending mode will result in underestimating N’ /M by 8 percent and C x by 
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12 percent. If the first four bending modes were included, then this 
error would be reduced to underestimating N'/M by 7 percent and C x by 
8 percent. 


Table II 

Comparisons of Aeroelastic Coefficients Resulting From A 
Polynomial Curve Fit Versus Bending Mode 



N' /M 

M/sec 2 

1 /rad 

n 

rad 

sec 2 

1 /rad 


Rigid 

Bending 

Total 

Aeroelastic 

Rigid 

Bend ing 

Total 

Aeroelastic 

Exact Solution 

6.72 

7.36 

7.52 

-0.174 

-0.233 

-0.247 

8 th Order Polynomial 

6.72 

7.36 

7.52 

-0.174 

-0.234 

-0.249 

1 st Bending Mode 

6.72 

6.92 

6.94 

-0.174 

-0.214 

-0.219 

1 st and 2 nd Bending 
Modes 

6.72 

6.94 

6.96 

-0.174 

-0.228 

-0.234 

ist s 2 nd, & 3 rd 
Bending Modes 

6.72 

6.97 

6.98 

-0.174 

-0.228 

-0.234 

1 st , 2 nd , 3 rd , & 4 th 
Bending Modes 

6.72 

7.00 

7.01 

-0.174 

-0.229 

-0.235 


Saturn V Flight Conditions; Mach =1.2 

a = 12° 

q = 30,000 N/M 2 


IV. CONCLUSIONS AND RECOMMENDATIONS 


The Saturn aeroelastic effects are significantly important during 
the period of maximum dynamic pressure. At this period the induced 
aerodynamic loading decreases the vehicle's stability and increases 
the component of local aerodynamic bending moment. The large external 
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body moments recorded as unexplained during Saturn I Block II flight 
tests result from changes in aerodynamic load distributions induced by 
aeroelasticity. 

It is recommended that six-degrees -of-freedom simulation programs, 
such as prelaunch wind monitorship, be modified to include the Saturn 
aeroelastic loadings. Such a program modification could be readily 
incorporated by storing the precomputed aeroelastic loads as a function 
of Mach number, dynamic pressure, and angle of attack. While this 
represents an appreciable amount of information, it would need to be 
computed only once for any geometric and stiffness configuration of the 
Saturn vehicle. An alternate approach would be to precompute the aero- 
elastic effects for specific trajectories and alter this input for 
significant changes in flight trajectories. 

An investigation of the downwash effects, which were neglected in 
this report, is recommended. This would require modifying the aero- 
dynamic loadings for the various elastic angle-of-attack distributions 
and using these loadings to determine the aeroelastic superposition 
solution. The investigation would be greatly enhanced by wind tunnel 
tests of Saturn models with known rigid deflected shapes. 
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FIG. 1. RIGID AND INDUCED 
AEROELASTIC NORMAL FORCE COEFFICIENT GRADIENTS 
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ACn a (1/deg) 


FIG. 2. INCREASE IN NORMAL FORCE COEFFICIENT 
GRADIENT DUE TO AEROELASTIC LOAD GROWTH 
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FIRST STAGE 
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SECOND STAGE 
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15° 00 
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FIG. 4. GEOMETRY OF 
THE APOLLO - SATURN FLIGHT VEHICLE 
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r IG. 5. INCREASE IN LOCAL AERODYNAMIC BENDING MOMENT DUE TO 
AEROELASTIC LOAD GROWTH FOR SATURN IB AT q * 20,000 N/m 2 




Increase in local Aerodynamic Bending Moment 
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AEROELASTIC LOAD GROWTH FOR SATURN IB AT q = 30,000 N/rn 









% Incease in Local Aerodynamic Bending Moment 



FIG. 7 EFFECT OF NONLINEAR ANGLE OF ATTACK 
AERODYNAMICS ON AEROELASTIC BENDING MOMENT 
AT STATION 1663 INCHES FOR MACH 1.0 
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% Increase In Local Aerodynamic Bending Moment 



0 10,000 20,000 30,000 40,000 50,000 


Dynamic Pressure (N/m 2 ) 
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FIG. 9. INCREASE IN LOCAL AERODYNAMIC BENDING 
MOMENT DUE TO AEROELASTIC LOAD GROWTH 
FOR SATURN Z AT STATION 3100 INCHES 
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Total Aeroelastic Effect (Bending 8 Load Growth) 

Bending Effect 

Rigid Body Aerodynamics 


C-i , Aerodynamic Moment Coefficient (rad/sec 2 )(1/rod) 



0 4 8 12 16 20 

Angle of Attack (deg) 


FIG. 12. SATURN I INCREASES IN 
AERODYNAMIC MOMENT COEFFICIENT AT MACH 1.2 
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Bending Effect 

Rigid Body Aerodynamics 


N'/m, Aerodynamic Force Coefficient (m/sec 2 ) (1/rad) 



FIG. 13. SATURN I INCREASES IN 
NORMAL FORCE COEFFICIENT AT MACH 1.7 
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FIG. 14. SATURN 1 INCREASES IN 
AERODYNAMIC MOMENT COEFFICIENT AT MACH 1.7 
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